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ABSTRACT 

It is of great interest to measure the properties of substructures in dark matter halos at galactic 
and cluster scales. Here we suggest a method to constrain substructure properties using the 
variance of weak gravitational flexion in a galaxy-galaxy lensing context. We show the effec- 
tiveness of flexion variance in measuring substructures in N-body simulations of dark matter 
halos, and present the expected galaxy-galaxy lensing signals. We show the insensitivity of 
the method to the overall galaxy halo mass, and predict the method's signal-to-noise for a 
space-based all-sky survey, showing that the presence of substructure down to 10 9 M© halos 
can be reliably detected. 

Key words: cosmology: observations - dark matter - gravitational lensing. 



1 INTRODUCTION 

The amount of substructure in dark matter halos on galactic and 
cluster scales is a question of considerable interest. ACDM cos- 
mological simulations predict thousands of dark satellites within 
the virial radius of the Milky Way, yet to date only ~ 20 have 
been observed (e.g. Kauffmann et al. 1993; Klypin et al. 1999; 
Diemand et al. 2007). Similar results have been obtained for our 
nearest neighbour Andromeda (McConnachie & Irwin 2006), and 
in galaxy groups (D'Onghia et al. 2007). Averaging over galaxies 
observed in the field, there appears to be a suppression in the ex- 
pected number below a baryonic mass of ~ 10 10 Mq (e.g. Read & 
Trentham 2005). This dearth of low mass galaxies could be telling 
us that galaxy formation becomes increasingly inefficient below the 
peak of the luminosity function (e.g. Dekel & Silk 1986; Bullock 
et al. 2000; Benson et al. 2002; Kravtsov et al. 2004; Read et al. 
2006). Alternatively, it could be telling us something about the na- 
ture of dark matter (e.g. Bode et al. 2001; Barkana et al. 2001; 
Avila-Reese et al. 2001; Knebe et al. 2008), or about the details of 
inflation (Zentner & Bullock 2003). 

Gravitational lensing is a powerful tool for probing substruc- 
ture within galaxies and clusters. Unlike photometric galaxy sur- 
veys, lensing directly probes the dark Universe; in principle, even 
if completely dark subhalos exist, these will have some influence on 
the observed gravitational lensing signal and may therefore be de- 
tected. Several lensing techniques for constraining substructure al- 
ready exist: one can examine the strong lensing of quasars by galax- 
ies and clusters (e.g. Mao & Schneider 1998; Metcalf & Madau 
2001; Metcalf & Zhao 2002; Keeton et al. 2003; Moustakas & 
Metcalf 2003; Kochanek & Dalai 2004; Metcalf et al. 2004; Met- 
calf 2005; Aazami & Natarajan 2006; Amara et al. 2006; Miranda 



& Maccio 2007; Shin & Evans 2008) including the time delay 
phenomenon (Keeton & Moustakas 2009). One can also constrain 
the substructure using weak galaxy-galaxy shear (Natarajan et al. 
2004). Here we examine the usefulness for measuring substructure 
of another weak lensing phenomenon: flexion, building on several 
earlier studies (c.f. Goldberg & Natarajan 2002; Goldberg & Bacon 
2005; Bacon et al. 2006; Irwin & Shmakova 2006; Schneider & Er 
2008). 

It has been noted by these authors that flexion responds to 
small-scale variations in the gravitational potential. In a galaxy- 
galaxy lensing context, the mean flexion in annuli around a halo 
will fall off rapidly, as the mean halo density gradient is small 
away from the central region. However, if substructure is present, 
the flexion variance in annuli may not fall off so quickly, as the 
substructures will lead to potential fluctuations which will cause a 
flexion varying rapidly from place to place. It is this idea that the 
current paper will examine. 

The topic is of particular interest in the context of forth- 
coming and planned large lensing surveys. Ground-based surveys 
such as Pan-STARRS 1 and the Dark Energy Survey 2 will obtain 
many thousands of square degrees of lensing-quality data, allowing 
very precise galaxy-galaxy lensing constraints. Beyond this, space- 
based survey telescopes such as Euclid 3 will obtain extraordinary 
accuracy for galaxy-galaxy lensing due to a high number density 
of galaxies and a survey area of 20000 square degrees. This pa- 



1 http://pan-starrs. ifa.hawaii.edu 

2 https://www.darkenergysurvey.org 

3 http://www.ias.u-psud.fr/imEuclid 
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per will make some initial predictions for the level of accuracy of 
flexion variance measurements with the latter survey. 

The paper is organised as follows. In section 2, we consider 
the relevant theory of weak gravitational lensing, concentrating on 
the flexion formalism. We proceed to define the galaxy-galaxy flex- 
ion statistics that we require, and show how they can be estimated 
theoretically. In section 3 we describe N-body simulations of halos 
with and without substructure, and explain the procedure for cre- 
ating lensing maps from these simulations. In section 4 we use the 
maps to illustrate the sensitivity of flexion to substructure, then cal- 
culate the expected galaxy-galaxy flexion variance signal together 
with its signal to noise for a Euclid-like survey. We examine the 
effect of varying the dominant halo mass, finding little impact, and 
show the much more pronounced effect of varying the amount of 
substructure present. This latter finding will illustrate the useful- 
ness of the technique as a probe of dark matter temperature. We 
present our conclusions in section 5. 



2 THEORY 
2.1 Flexion 

The flexion formalism described here is more fully developed in 
Bacon et al. (2006). We begin by noting that lensing in the weak 
regime can be described by mapping the surface brightness of a 
galaxy in the source plane, fs(/3i), to the surface brightness in the 
image plane, fi(6i): 



fi(0i) = fs(Pi) = fs Airfi + -I );„'), Ik 
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(1) 



Here we have introduced several quantities: firstly, the Jacobian 
matrix A, which is taken to be constant across the galaxy image 
in the weak regime. It can be written in terms of lensing quantities 
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where k is the convergence; this maps a unit circle in the image 
plane to a circle with radius 1 — « in the source plane. In the case 
of an isolated lens, the convergence is proportional to the projected 
density of matter in the lens (c.f. Bartelmann & Schneider 2001), 



where E is the 2D projection of the density p, 
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where the integration is over the radial distance 7-3 and E cr is given 
by 
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where Di, D s and Di s are the angular diameter distances from the 
observer to the lens, from the observer to the source, and from the 
lens to the source respectively. If we also define the lensing poten- 
tial ip, which is proportional to the projection of the gravitational 
potential $, 

4vrG DiD, 
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then we can also write the convergence as 
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The other term in A is the shear -y,; this maps a circle in the im- 
age plane to an ellipse in the source plane. Its components can be 
written as 
1 

71 = 2 

The next term in equation (1) is the 77-tensor; this contains the 
lensing information at the next order of approximation, and corre- 
sponds to the varying of convergence and shear across an object. 
As shown in Bacon et al. (2006), the D-tensor can be written in 
terms of the flexions, 
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where Gi are the components of 3-flexion, describing the degree to 
which an object resembles a trefoil, and Fi are the components of 
1-flexion, describing the skewed shape of an object. We will only 
consider the 1-flexion for the purposes of this paper, which has the 
property of being the gradient of the convergence, 
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It is this property of 1-flexion that is so important for the technique 
of this paper; 1-flexion will respond wherever the density is varying 
rapidly from place to place, which is the case with substructure. 

2.2 Galaxy-Galaxy Flexion 

We can now introduce the flexion variance in annuli as a substruc- 
ture probe. It is usual in galaxy-galaxy lensing (e.g. Kleinheinrich 
et al. 2006; Mandelbaum et al. 2006) to choose particular galaxy 
samples which act as lens and source sets. The sets may have some 
members in common, depending on the technique used for lens- 
source correlation, but at any rate the selections can be achieved 
either via photometric redshifts, spectroscopic redshifts or a com- 
bination of the two. 

Let us orientate ourselves with the most familiar case of 
galaxy-galaxy lensing, which involves shears. The angular separa- 
tion 9 of foreground object / and background object b is measured, 
and the shear of b is decomposed into components tangential and 
diagonal to the line connecting / and b, 

jt = -» ( 7 e~ 2 ^) IB = [ie- 21 ^ (1 1) 

where (j> is the position angle of the line, 7 = 71 + 172, and 5R 
and 9 denote real and imaginary parts respectively. For a circular 
foreground lens, jt will be activated by gravitational lensing, while 
7_s will not; if it is present, it is due to systematic effects. 

Similarly a background flexion can be decomposed, but 
whereas shear is decomposed into tangential and diagonal compo- 
nents, flexion is decomposed into radial and tangential components 
due to the fact that it has different rotational properties to shear: 
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where F = F\ + iF 2 . Here the r component is activated by gravity 
for a circular lens, while the B component is unactivated unless 
systematics are present. 

The behaviour of flexion at an angular distance 9 from a singu- 
lar isothermal sphere (SIS) was considered by Bacon et al. (2006). 
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They show that while the surface density of the SIS is proportional 
to O^ 1 , the 1-flexion and 3-flexion drop as 9~ 2 . While real galaxies 
are not truly SIS, they are close enough to this profile that if we take 
foreground-background pairs, measure F r , and consider the mean 
F r in annuli, we will see a similarly rapid drop in galaxy-galaxy 
flexion with annulus radius 8. 

This will remain true if there are substructures at angular dis- 
tance 6 from the centre of foreground galaxies. Although in the 
extreme locality there will be a larger flexion than usual, the mean 
signal averaged around the galaxy will still drop rapidly with 9. 

However, this will not be true for the flexion variance in the 
annulus. This will respond to any density fluctuations within the 
annulus. Therefore, in any annuli with non-negligible substructure, 
even if the mean flexion is small (as the mean gradient of density 
is small), the flexion variance will remain comparatively large. It is 
this behaviour that we will use to constrain substructure on galactic 
scales. 

At this point one may ask whether flexion variance is the best 
tool for our task; wouldn't flexion correlation functions in annuli 
provide more information? The question could be informed by ex- 
perience in cosmic shear studies, where shear correlation functions 
provide more finesse than shear variance in cells. 

However, in our present case a correlation function does not 
seem to be helpful. Since the correlation function in question would 
be in annuli around a foreground galaxy, it constitutes a form of 
galaxy-galaxy-galaxy lensing (Schneider & Watts 2005); it is a 
three-point statistic. In order to use this to measure substructures, 
it is necessary for two background sources to be close to the same 
substructure as well as to the foreground galaxy; this rarely hap- 
pens, leading to low signal-to-noise. On the other hand, the flex- 
ion variance in annuli only involves two points, a foreground- 
background pair, with much greater signal-to-noise as we shall see. 



3 SIMULATIONS 

In this section we describe N-body simulations which we will use 
to demonstrate the utility of flexion variance as a probe of substruc- 
ture. 

3.1 3-D Density 

We use the cosmological ACDM simulation already presented in 
Diemand et al. (2005). The simulation was run using PkdGRAV 
(Stadel 2001), with cosmological parameters: (Jl m , Qa, <Ts,h) = 
(0.268,0.732,0.7,0.71), and a box of size i box = 90 Mpc, 
with 300 3 particles. The initial conditions were generated with 
GRAFIC2 (Bertschinger 2001). From the simulation volume, we 
extracted four Milky Way sized halos at a mass resolution of m p = 
5.7 x 10 5 M o ; their virial masses are [2.1, 1.5, 1.2, 1.3] x 10 12 M Q . 
While we are therefore very limited in our number of lenses (to 
three projections of each of four high-resolution galaxies), we will 
find that this is sufficient to give the initial indicative results re- 
quired by this paper. 

As in Read et al. (2008), the subhalos inside each 'Milky Way' 
at redshift z — were identified using the AHF algorithm (Gill 
et al. 2004). We considered all subhalos with > 50 particles and as- 
signed particles to the smallest structure they appear in so that each 
particle was counted only once. In some cases we will remove sub- 
structure; this is achieved by subtracting all particles not assigned 
to the main halo. An example halo, with and without substructure, 
is shown in Figure 1 . 




-200100 100 200 
x(kpc) 




-200-100 100 200 
x(kpc) 

Figure 1. N-body simulation of one of our Milky Way mass halos, 
with substructure (top panel) and without substructure (bottom panel), 
with contours showing projected logarithmic density over the range 
[0.02,672] M pc- 2 . 

3.2 2-D Convergence 

The 3D numerical simulations discussed above represent the den- 
sity field using discrete particles. We transform these into conver- 
gence maps by projecting the particles along particular spatial di- 
rections and placing the particles onto a 2D 1024x1024 grid, which 
we carry out using the IDL cloud-in-cell routine available as part of 
The IDL Astronomy User's Library 4 . We produce three projection 
maps for each 3D halo (by projecting along the x, y or z axis). 

We investigated a number of techniques for filtering the mass 
maps. This is important because the finite number of simulation 
particles introduces shot noise into the 2D maps; this can compete 

4 http://idlastro.gsfc.nasa.gov 
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Figure 2. Standard deviation of 1 -flexion in annuli for one galaxy halo, 
in units of (arcsec) -1 . Arrows show the positions of subhalos, with arrow 
length proportional to subhalo mass. 

with the substructure signal we are investigating. Here we show the 
results obtained using Multiscale Entropy Filtering (MEF) (Starck 
et al. 2006). This provides superior performance to a simple Guas- 
sian filter, as it reduces the shot noise while preserving the den- 
sity fluctuations, as we shall demonstrate. For this purpose we use 
the routines supplied in the public code MRLens (Multi-Resolution 
methods for gravitational LENSing) 5 , choosing the eight scale for 
MEF with first detection at the second scale. The 2D mass maps (E) 
are then rescaled using the critical density (E c ) to give the conver- 
gence (k). We will use a lens-source configuration of D s = 1200 
Mpc, Di — 860 Mpc, and Di s = 610 Mpc corresponding to our 
choice of median redshifts for foreground (z = 0.46) and back- 
ground (z — 1.1) sources appropriate for our projected Euclid sur- 
vey. 



4 ANALYSIS 

Armed with this set of simulations, we are in a position to examine 
whether substructure can be reliably measured using galaxy-galaxy 
flexion variance. 

4.1 Flexion Variance Sensitivity to Substructure 

Using the convergence maps for each halo, we calculate the related 
1 -flexion map using Fi = din, and smooth the resulting flexion 
map with a 0.6" diameter top hat filter to remove small-scale flex- 
ion peaks inaccessible to galaxy shape measurements. We then cal- 
culate the mean and standard deviation of F in annuli with width 
0.8", centred on the mode of the galaxy's k distribution. 

As an initial example, Figure 2 shows the standard deviation 
of for the radial 1-flexion in annuli, for one example halo. Also 
displayed are arrows showing the radial positions of subhalos found 

5 http : //irfu. cea.fr/Phocea/Vie_des_labos/Ast/astvisu.php7id_ast 
878 
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Figure 3. Mean (dotted lines) and standard deviation (solid lines) of conver- 
gence k in annuli for the twelve halo orientations combined. Upper lines in 
each case show results with substructure, while lower lines show the results 
when substructure is removed. 



in the simulation, with arrow size proportional to the mass of the 
subhalo. We can see that there is a rather close correspondence be- 
tween the flexion variance and the subhalo positions and masses. 
This encourages us to examine what the signal will be for an en- 
semble of galaxy halos in a galaxy-galaxy lensing context. 



4.2 Galaxy-Galaxy Flexion Signal 

In this section we will present an estimate of the galaxy-galaxy 
flexion signal; but in order to begin with a familiar quantity, we 
first examine the convergence for our halos in annuli of width 0.8", 
centred on the mode of the main halos' n distribution. Figure 3 
shows the mean convergence R and standard deviation <j k for an 
ensemble of background galaxies behind the twelve stacked halos, 
in annuli with radius 9. Note that a K is smaller than R for 9 < 
50"; since k is proportional to the projected surface density, this 
reflects the fact that the substructure fluctations are fairly small in 
amplitude in relation to the mean density. Hence the difficulty in 
using convergence or shear to measure substructures; the signal-to- 
noise on a K will be smaller than that on R. 

The upper and lower lines for R and a K can be compared to 
see the effect of including or omitting substructure; the mean of the 
convergence is a little higher on scales > 30" if we include sub- 
structure, and the standard deviation of the convergence is several 
times larger. This is to be expected; the presence of substructure 
changes the profile a little, and significantly alters the spatial varia- 
tion of the matter distribution. 

We now consider the 1-flexion galaxy-galaxy signals, again 
calculated in annuli with width 0.8", centred on the mode of the 
galaxy's k distribution, and smoothed with a 0.6" diameter top hat. 
Figure 4 shows the mean F and standard deviation of for the radial 
1-flexion in annuli, for an ensemble of background galaxies behind 
the halos. We show results with either MEF filtering or no filtering, 
and with or without substructure, as it is important to understand 
the impact of our filter on the results. We should emphasize that this 
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Figure 4. Mean (dotted lines) and standard deviation (solid lines) of 1- 
flexion in annuli for the twelve halo orientations combined, in units of 
(arcsec) -1 . Upper lines in each case show results with substructure, while 
lower lines show results when substructure is removed. Top panel: MEF 
results; bottom panel: unfiltered results. 



filter is applied to the N-body simulations to increase their realism; 
it is not a filter applied at the measurement stage. 

We immediately see that the filtering of the N-body simulation 
is valuable; in the cases where no filtering is applied, the flexion 
standard deviation is substantial even in the absence of substruc- 
ture. This is due to the particle noise in the simulations that leads to 
spurious high frequency convergence gradients. However, with the 
MEF filter in the absence of substructure, the flexion standard de- 
viation is at least an order of magnitude smaller than the unfiltered 
standard deviation. 

Importantly, we see that when substructure is included (with 
MEF filtering of the simulations), o dominates over F on scales 
greater than about 10", and is two orders of magnitude larger than 
the no-substructure signal. The former observation confirms our 
claim that flexion variance is of interest, as it has much larger 



Figure 5. Ratio of oj, with modified central halo amplitude to the usual 
Opt, Solid: no central halo; dashed: total mass standard deviation of 4 X 
10 12 A/q; dotted: total mass standard deviation of 1O 13 M0. 

signal-to-noise than standard galaxy-galaxy flexion. The latter ob- 
servation shows that it is substructure which is activating this sig- 
nal; this is important, as the overall ellipticity of the central halo 
might have given a flexion variance (c.f. Hawken & Bridle 2009) 
but our result demonstrates that this is negligible for our statistic. 
To measure halo ellipticity with flexion, one should instead follow 
the methods of Hawken & Bridle (2009). 

The flexion variance also dies off more slowly with than 
F does. This is presumably due to the different phenomena being 
probed: F is probing the mean gradient of the overall galaxy halo 
density; while op is dominated by the gradients of substructure 
halos within the annuli, which remain similar regardless of which 
annulus they are found in. 

4.3 Insensitivity to Varying Dominant Halo Mass 

An issue that we have not yet dealt with is the question of whether 
combining galaxies with various masses contributes a dominant un- 
wanted signal to our galaxy-galaxy flexion signal. It is common in 
galaxy-galaxy lensing to divide up or scale one's foreground lens 
set according to a mass proxy such as luminosity (e.g. Kleinhein- 
rich et al. 2006) or stellar mass (e.g. Mandelbaum et al. 2006). Nev- 
ertheless, within such a bin in mass, will the flexion variance signal 
be dominated by the overall mass variance rather than the substruc- 
ture signal? 

In order to test this, we measured the mean flexion signal for 
each of the main halos. We then subtracted this signal from each 
halo's flexion map, leaving the flexion due to substructure alone. 
We then optionally re-added a random proportion of the mean flex- 
ion signal, leading to a total mass standard deviation of 4 x 10 12 M@ 
or 1O 13 M0, and remeasured the ensemble flexion variance signal 
for our 12 realizations. The effects on our signal are shown in figure 
5. 

Firstly we notice that the difference between the shear vari- 
ance with and without central halos is very small; most of the sig- 
nal in our original ensemble is therefore coming from substructure. 
This is borne out by the case where we have a halo mass variance 
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of 4 x 10 Mq. The ratio of this signal to our original is within a 
few per cent of unity beyond 8 > 20" . In the extreme case where 
the mass variance is 1O 13 M , we obtain about 10% of the flexion 
variance from the main halosforfl > 30 , but this case corresponds 
to hardly having any mass binning at all. 

We conclude, then, that provided reasonable steps are taken 
to deal with the galaxy-galaxy lensing in mass or luminosity bins, 
flexion variance is totally dominated by the substructure signal. 

4.4 Signal-to-Noise for Galaxy-Galaxy Flexion 

The question now arises whether flexion variance is measurable 
with sufficient accuracy on the relevant scales to constrain sub- 
structure. The difficulty is that source galaxies are measured with 
a substantial intrinsic flexion (c.f. Goldberg & Bacon 2005); the 
distribution has a strong central peak and wide wings, such that the 
measured distribution has a 68% range of — 0.1 arcsec -1 . Clearly 
this intrinsic variance dominates over the substructure variance at 
all scales. However, following a statistical approach from Bacon 
et al. (2000), one can estimate the flexion variance due to substruc- 
ture, 0"s„;,, as an excess variance within an annulus, over and above 
that due to shape noise cr 2 ntrinsic : 

® sub — ^annulus ^intrinsic (13) 

where cr^„„„ Ius is the measured total variance within an annulus, 
while ^Intrinsic ls the variance measured for the whole ensem- 
ble of source galaxies. Since this equation is only exact for Gaus- 
sian distributions, we have checked that it is approximately true for 
our simulation flexion distribution and realistic shape noise distri- 
bution; we find that a 2 ub is correctly estimated to within 5%. As 
described by Bacon et al. (2000), the error on this estimator is ap- 
proximately 

~ o- 2 ntrinsic /\fN (14) 

where N is the number of objects in the annulus. 

We can now use this equation to estimate how strong a signal 
we expect for a Euclid-like survey. We use OFint = 0.1, and calcu- 
late N using the standard survey parameters of Amara & Refregier 
(2007). We examine the case of choosing 9 foreground galaxies 
per sq arcmin with z me d = 0.46, 26 background galaxies per sq 
arcmin with z me <i = 1-1. and a survey area of 20000 sq deg. Fig- 
ure 6 shows the expected value of our estimator o Fsub (solid line) 
and the noise level on this estimator in each annulus (dashed line). 
Note that for a large range in 9, the S/N is considerable: S/N ~ 1 
to 5 for each of 50 annuli between 20" and 60". This suggests that 
substructure can be studied in great detail using this method. 

4.5 Sensitivity to Substructure Content 

To pursue this point, we provide some examples of the degree to 
which we can distinguish between different substructure scenarios. 

First, we examine the impact of removing half of the subhalos 
at random in each of our simulations. We carry out the analysis of 
sections 4.2 and 4.4 for these modified simulations, and show how 
their flexion variance signal differs from the usual case on Figure 
7. Here we have used 50 bins between and 80" for clarity. 

We note that the flexion variance is itself approximately 
halved in this case. A \ 2 measure for the significance of this differ- 
ence is the sum of (c F — a Fhal f) 2 /a[al ub ] 2 for all annuli, and is in 
this case ~ 250, confirming that such a substructure configuration 
would be strongly distinguished. 




theta (arcsec) 

Figure 6. Flexion variance estimator (solid line) in units of (arcsec) -2 , and 
expected noise level in each annulus (dashed line) for the Euclid survey 
described in the text. 
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Figure 7. Ratio of flexion variances with substructure and without 50% of 
the subhalos. Grey shading shows region within la uncertainty of flexion 
variance with substructure. 



However, that scenario is not expected physically. More plau- 
sibly, if dark matter has a non-negligible temperature, this will pref- 
erentially remove substructures up to a certain mass threshold. To 
examine whether this type of phenomenon would be detectable by 
flexion variance, we remove all halos with mass < 5 x 1O S M0 or 
< 10 9 A/ o from our simulations. We again carry out the analysis 
of sections 4.2 and 4.4 for the modified simulations. The results of 
this process are shown on Figure 8. 

We see that there is only a small difference between the sig- 
nals for all subhalos and for all subhalos > 5 x 10 8 Mq; they are 
only slightly distinguishable within the error expected as shown by 
grey shading, with X 2 = 4.3. Care should be taken in drawing con- 
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Figure 8. Ratio between (a) flexion variances with all subhalos, and with 
all subhalos > 5 X 10 8 Mq (dashed line); (b) flexion variances with all 
subhalos, and with all subhalos > 1O 9 M0 (solid line). Grey shading shows 
region within lcr uncertainty of flexion variance with all subhalos. 

elusions from this, as halos with mass 5 x 10 Mq are still near 
our resolution limit; we will explore lower mass substructure in fu- 
ture work. On the other hand, the difference between the > 10 9 M Q 
case and the full substructure case is easily detected, with \ 2 = 37. 
According to the thresholds given by Barkana et al. (2001), this 
would approximate detection of a 2keV mass scale for warm dark 
matter. 



5 CONCLUSIONS 

In this paper we have explored the utility of galaxy-galaxy flexion 
variance for the purpose of measuring the degree of substructure on 
galactic scales. 

We have described the relevant statistics, showing how flex- 
ion radial and systematic modes are constructed in a galaxy-galaxy 
lensing context, and introducing the concept of flexion variance in 
annuli. We have explained how this is a more suitable probe than 
galaxy-galaxy-galaxy flexion (or flexion correlation functions in 
annuli), as the latter has a very rapid drop in signal as a function 
of angle. 

We have gone on to test the use of these statistics by adopting 
a set of N-body simulations of galaxy halos, including substructure, 
or removing this substructure by means of halo-finding algorithms. 
We have then calculated the weak lensing convergence associated 
with these halos, applying the suitable MEF filter to reduce the ef- 
fect of particle shot noise. 

We have shown that the flexion variance is able to detect sub- 
structures on a halo by halo basis. This carries through to the full 
galaxy-galaxy flexion variance expected with an ensemble of ha- 
los; the signal found is substantially larger than that for the flexion 
mean, or indeed the flexion variance in the absence of substruc- 
tures. The underlying central halo mass variance is not found to be 
a dominant source of noise for this signal. 

We have made predictions for the level of the flexion vari- 
ance signal-to-noise for a space-based survey such as Euclid, find- 



ing that substructure amplitudes will be measured with significant 
precision. This allows us to discriminate between substructure sce- 
narios with different numbers and masses of halos, and will enable 
a constraint on dark matter particle mass. 
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